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A powerful new impurity solver is shown to permit a systematic study of the doping driven Mott 
transition in a one-band Hubbard model within the framework of single-site dynamical mean field 
theory. At small dopings and large interaction strengths we are able access low enough temperatures 
that a reliable extrapolation to temperature T — may be performed, and ground state energies of 
insulating and metallic states may be compared. We find that the T — doping-driven transition is 
of second order and is characterized by an interaction-strength dependent electronic compressibility, 
which vanishes at the critical interaction strength of the half filled model. Over wide parameter 
ranges the compressibility is substantially reduced relative to the non-interacting system. The metal 
insulator transition is characterized by the appearance of in-gap states, but these are relevant only 
for very low dopings of less than 3%. 



I. INTRODUCTION 

The 'Mott' or correlation-induced insulating state is 
a fundamental unifying concept in modern condensed 
matter physics. The physical properties of many inter- 
esting materials, including organic conductors^ colossal 
magnetoresistance manganites, 2 actinides such as Ce and 
Pu)&^ and many transition metal oxide compounds^ are 
believed to be controlled by proximity to a Mott insu- 
lating state. However, the accurate theoretical descrip- 
tion of the physics of Mott insulators poses challeng- 
ing problems, and many questions remain unresolved. 
A particularly important class of open questions, cru- 
cial, for example to the physics of high temperature 
superconductivity** concerns the behavior at strong in- 
teractions as the carrier concentration is varied away 
from the commensurate values at which Mott insulating 
behavior occurs. 

An important theoretical step forward was achieved 
with the development of "dynamical mean field theory" j£ 
which showed that if the momentum (p) dependence 
of the electron self energy S(p, oS) can be appropriately 
approximated, the computation of electronic properties 
may be reduced to the solution of a quantum impu- 
rity model, along with a self consistency condition. In 
"single-site" dynamical mean field theory, the momentum 
dependence is completely neglected, £(p, w) — > 
and the impurity model is a single site coupled to a 
free fermion bath. This approximation is strictly valid 
in the limit of inifinite coordination number, but cap- 
tures many features of the behavior of finite dimensional 
compounds fi^i^ 

One of the early successes of the single-site dynamical 
mean field theory was an explication of the basic phase 
diagram of the Mott transition in single band (Hubbard- 
like) models^ If magnetism may be neglected, the phase 
diagram in the space spanned by temperature (T), inter- 
action strength (U) and chemical potential (/z) takes the 
form shown in the upper panel of Fig. [TJ At low, but 
non-zero temperature, a first order metal-insulator tran- 
sition occurs on the surface delimited by thick red lines. 
The coexistence region, where both metallic and insulat- 
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FIG. 1: (color online) Sketch of the phase diagram of the 
Hubbard model in the single site DMFT approximation (semi- 
circular density of states of bandwidth 4t) with magnetic or- 
der suppressed by averaging over spin. Top panel: thick red 
lines indicate the surface in the space of temperature, interac- 
tion strength and chemical potential, where a first order metal 
insulator transition occurs. Thin black lines delimit the coex- 
istence region (hashed for cuts across the n = /ii and U — 6t 
planes) where both metallic and insulating solutions to the 
DMFT equations exist. Lower panel: T — phase diagram 
and spectroscopic gap A. 



ing solutions to the DMFT equations exist, is indicated 
by hashed areas on the planes /i = fii(U) (corresponding 
to half-filling) and U = 6t. The T = phase diagram 
is shown in the lower panel of Fig. [TJ If the chemical 
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potential is held at the value /j,x(U), then the model is 
metallic for interactions weaker than a critical value con- 
ventionally denoted U C 2 and we expect properties vary- 
ing smoothly with chemical potential. As U — > U C 2 from 
below, the quasiparticle weight and compressibility van- 
ish continuously, so that the metal-insulator transition 
at T = is second order. For U > U C 2 the model is 
insulating for carrier concentration n — 0.5 per spin, and 
there is correspondingly a region of chemical potential 
(bounded by the curve labeled fi C 2)> where the density is 
pinned at 2n = 1 and no metallic solution exists. 

The second order character of the transition at U C 2 
strongly suggests that the critical chemical potential [i c2 
smoothly approaches fix as U —> U C 2 from above. How- 
ever, the insulating phase at U > U C 2 is characterized by 
a spectroscopic gap A which does not vanish at U = 
and remains locally stable for a range of U < U C 2 (the 
lower limit of the region in which the insulating solution 
exists, is conventionally denoted by U c i). The gap, pre- 
sented as a chemical potential difference from the half 
filled value fix is shown as a dashed-dotted line in Fig. [1] 
The difference between the values of A and \x C 2 implies 
that in this approximation, doping a Mott insulator pro- 
duces "in-gap" states, and Fisher, Kotliar and Moeller 
have presented more precise analytical arguments which 
support these ideas and show that doping generates in- 
gap states for all U > U C 2&^ 

However, moving beyond general arguments, very little 
is known with confidence about the specifics of the phase 
diagram. The essential difficulty has been the lack of 
numerical methods powerful enough to address the region 
of strong correlations, low temperature and low doping 
where the interesting physics occurs. In this paper we 
use a newly developed metho d 11 ' 12 to solve the problem. 
The method permits access to strong correlations and 
low temperatures, with an unprecedented accuracy which 
enables us to construct thermodynamic potential curves 
and establish the nature and location of the transition. 
Our results are consistent with the following scenario: the 
point U = U C 2 and fi = fix is a quantum critical point 
at which the electronic compressibility (proportional to 
dn/dfi) vanishes linearly in \fi — fi\\. For U > U C 2 and 
T = 0, a second order metal-insulator transition occurs 
at an interaction dependent chemical potential \i C 2 (U) — 
fix ~ (U — U C 2) X , with an exponent x close to 1/2. For 
U > U C 2 the compressibility does not vanish as /i — > fi C 2, 
but at large enough fi — fix, dn/dfi ~ | /x — A*x | - The 
physics associated with the critical point is visible over a 
reasonable range of U > U C 2- 



II. MODEL AND FORMALISM 
A. Model 

In this paper we present results for the paradigm 
strongly correlated model, the one-orbital Hubbard 



model, defined on a lattice of sites i by 

h = - ^kW^- + Yl Un ^ n i,i- (!) 

i,S,a i 

The energy dispersion e p is defined as the Fourier trans- 
form of t(5) and the only property of the dispersion 
which will be important for us is the density of states 
V(lo) — J^p^i 10 ~ £ p)- I 11 our specific calculations we 
shall take 

V(lo) = V4i 2 -w 2 /(27rf 2 ). (2) 

For this choice of density of states the chemical po- 
tential /i corresponding to the (potentially Mott insulat- 
ing) density of one electron per site is fix = U/2 and 
the critical interaction strength for the zero temperature 
Mott transition in the single-site DMFT approximation 
is U c2 w 5.8*^ 

B. DMFT method 

The single-site DMFT reduces the solution of the lat- 
tice problem to the solution of a quantum impurity prob- 
lem defined by 

Hqi = -^2ndld a + Un dA n dd 

a 

+ J2 (VedUe,* + V e *cl a d a + ecl a c^) . (3) 

It is useful to define the hybridization function^ 

Fa(r) - l^l 2 (T r ct iff (r)c e , CT (0)) bath , (4) 

e 

the d-electron Green function G a (r) = — (T T d rr (T)d a (0)) 
and self energy E CT = d T + \i + F a (j3 — r) — G^t) -1 . 
The hybridization function is fixed by the self consistency 
condition 

G a (u) = [ de V{£) (5) 

J u - e- n-Y, a \u}) 

The challenging numerical task is computing G(t). 
In order to do this we have recently developed a new 
solver , 11 ' 12 which is based on a diagrammatic expansion 
of the partition function in the impurity-bath hybridiza- 
tions and the Monte Carlo sampling of certain collec- 
tions of the resulting diagrams. The summation of di- 
agrams into determinants eliminates the sign problem, 
even away from half-filling, and our approach, which ex- 
pands around an exactly solved atomic limit, leads to 
lower perturbation orders at stronger interactions U. The 
method thus allows unprecedented access to low temper- 
atures and strong interactions^ and will be used here to 
study the doping-dependent Mott transition. Near the 
end points r = and /3 the Green function converges 
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very rapidly, but more effort is needed to accurately de- 
termine the long-time behavior. We chose a resolution 
of 10000 points for the Green function and a smoothing 
procedure (averaging over 30 neighboring bins) at inter- 
mediate r to reduce the statistical errors in the region 
where such a high resolution is not necessary. With this 
resolution, the systematic errors should be small and we 
therefore estimated the error bars on quantities such as 
densities and energies from their variation in successive 
iterations of the converged solution. Where no error bars 
are given, the errors are smaller than the symbol size. 

From the computed G(r) we directly obtain the density 
per spin n = G(r — ► while the "internal energy" may 
be computed as 

E(fx,T) = (H-nN) 

= 2t 2 / (ItG(t)G(-t) +UD- 2(in. (6) 
Jo 

Here we have used a property of the semicircular density 
of states to obtain a compact expression for the kinetic 
energy term, while the expectation value of the interac- 
tion term is obtained from a direct measurement of the 
double occupancy D. 

C. Extrapolation to T = 

As will be shown below, a characterization of the 
metal-insulator transition requires the construction of 
thermodynamic potentials and the extrapolation of our 
data to T = 0. The insulating state is characterized by a 
large gap (which means that the energy at low tempera- 
ture is exponentially close to the ground state value) and 
an extensive spin entropy of In 2 per site, so the thermo- 
dynamic potential of the insulating state is 

n ins = E inB -T\n(2), (7) 

where E lns is computed from Eq. ©. 

The entropy of the metallic state is in general not easy 
to obtain. However we note that within single-site dy- 
namical mean field theory the metallic phase is, at low 
temperatures, a Fermi liquid characterized by a T 2 vari- 
ation of physical quantities. In particular, at sufficiently 
low T, the energy and thermodynamic potential of the 
metallic state are given respectively by 

E mct {n, T) = E met {n, T = 0) + i 7 T 2 , (8) 

IU0», T) = E met (», T = 0) — i 7 T 2 , (9) 

where the specific heat coefficient 7 = limx-^o C/T is 

( 10 ) 

Here /xq is the chemical potential which, in the model 
with U = 0, produces the density corresponding to /1 



in the interacting model. In the last line we used the 
fact that for the dopings of interest the density of states 
may be approximated by its half-filled value and that the 
self energy derivative at low T may be approximated by 
the value of the imaginary part of the self energy at the 
lowest Matsubara frequency divided by -kT . 

We estimate 7 in two ways: from the electron self en- 
ergy, via Eq. (flO]) , or by fitting the measured energies to 
a T 2 dependence and using Eq. (jHJ). For \i — » ^ c2 (espe- 
cially near XJ c i) the range over which the energy obeys a 
T 2 law becomes small and we find that obtaining 7 from 
the self energy leads to smaller errors, which we estimate 
to be at the 10% level. We have verified that for the dop- 
ings considered, we can reach low enough temperatures 
that the metallic entropy S = 7T is much smaller than 
In 2. 

Finally, we note that the thermodynamic potential 
may alternatively be obtained from the density-chemical 
potential trace via the thermodynamic relation 2n — 
—dCl/d/i. Choosing a reference chemical potential \i c we 
have 

fimet(M) = ^mot(Mc) - 2 ( dn'n(n'). (11) 

J He 

Below, we will use the T = limit of Eq. (fill) to 
demonstrate the consistency of our analysis and show 
that at T — the doping driven transition becomes sec- 
ond order. 



III. RESULTS 
A. First order transition at T > 

The upper panel of Fig. [2] shows the variation with 
chemical potential, fi, of the carrier concentration per 
spin, 71, measured relative to the Mott insulating value 
0.5 at the very low temperature t/T = j3t = 400. For the 
weakest interaction strength, U = 5.6t, the carrier con- 
centration varies smoothly with chemical potential, im- 
plying that the phase is metallic even at half filling. For 
the larger interaction strengths a gap (region where n is 
approximately independent of /i) is visible, showing that 
for these [/-values the model at (3t — 400 is insulating 
for a range of chemical potentials. At this temperature 
the critical interaction strength for the Mott transition is 
therefore U/t~ 5.65, consistent with the accepted phase 
diagram of the half filled models 

The density-chemical potential traces shown in the 
lower panel of Fig. [2] highlight an unusual scaling be- 
havior near the critical point at U = U C 2 and n = 0.5 
(ji = Hi = U/2): the density (measured from 0.5) varies 
as the square of the chemical potential (measured from 
/ii), in other words, near the Mott point, the compress- 
ibility per spin, which is up to a factor l/(2n) 2 given 
by dn/d/i, vanishes proportionally to |/i — with a co- 
efficient 0.44//ii 2 . The U — 5.6i curve exhibits at the 
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FIG. 2: Doping per spin, 0.5 — n, as a function of chemical 
potential for j3t = 400 and indicated values of U/t. At this 
temperature, the transition at half-filling (/z = fix = U/2) 
occurs at U C (T) fn 5.65. For larger interactions, a gap opens 
and shifting the chemical potential induces a first order metal- 
insulator transition. Near the critical point U C 2, n — 0.5 ~ 
(fi — /ii) 2 , but as one moves away from the critical point, the 
onset of doping becomes linear in /i. 
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FIG. 3: Doping- vs-ji for U/t = 6 (upper panel), 6.5 (lower 
panel) and fit = 100, 200 and 400. For fixed fx, the dop- 
ing increases with decreasing temperature. The vertical lines 
indicate the positions of the first order phase transition, as 
estimated from the crossing of the thermodynamic potential 
curves (dashed: fit = 200, solid: fit = 400). Solid trian- 
gles show the densities obtained from the thermodynamic po- 
tential using 2n — —dQ,/dfi (up-triangles: fit = 200, down- 
triangles: fit = 400). 



smallest [i a crossover away from the square root behav- 
ior to the constant dn/dfj, expected in a metallic phase. 

For U larger than the critical value, the curves exhibit 
a slight downward trend away from the (/x — fXi) 2 scaling, 
indicating a linear onset at very small dopings, but more 
importantly the curves are cut off by a discontinuity in- 
dicating our inability to numerically stabilize a metallic 
phase, and suggesting that at T > the doping driven 
metal-insulator transition is first order. The first order 
transition also occurs in the interaction driven (n = 0.5) 
case£ and is expected from the extensive entropy (In 2 
per site) of the paramagnetic insulating state. 

We now consider in more detail the behavior at very 
small dopings. The two panels of Fig. [3] show density- 
chemical potential traces for three low temperatures 
T/t = 0.01, 0.005 and 0.0025. These results show that for 
dopings (per spin) > 0.01, the lowest temperature data 
are essentially converged to the T = result, while at the 
smallest doping some temperature dependence clearly re- 



mains. 

The data presented in Fig. [3] define the range of pa- 
rameters over which the metallic state can be stabilized 
by our numerical procedure. The jump in density sug- 
gests the presence, in the T > phase diagram, of a 
first order metal-insulator transition, but the computed 
jump position is a spinodal point. To locate the first 
order phase transition, we compute the thermodynamic 
potentials f2i, ls and r2 mct using Eqs. and ^ and 
the specific heat coefficients obtained from the self en- 
ergies. Representative examples are shown in Fig. [4] 
which plots an approximation to the quasiparticle weight 
Z = 1/(1 - dZ/duj\^o) « 1/(1 - 3mE(w )/wo) as a 
function of doping for U/t = 6,6.5 and 8 and several 
chemical potentials. Note that if we define total doping 
x = 2(0.5 — n), then our data for U/t — 8 are roughly 
consistent with jt = 1.9/x, those at U/t — 6.5 with 
-ft = 1.35/x and those at U/t = 6 with -ft — 1.0/x. 
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FIG. 4: Quasi-particle weights estimated using Z « 1/(1 — 
S$mE(cJo) / 'wo) , evaluated at the lowest Matsubara frequency 
u>o, plotted as a function of doping per spin. The data points 
connected by dashed lines correspond to /3t = 400 and ap- 
proximation (|10[) . Solid dots show an estimate for Z which is 
based on a 2-parameter fitting function. 



Especially near U c %, the measured Z do not quite ex- 
trapolate to zero as i — > 0, which may be due to the 
approximation of the derivative in Eq. (I10p . The full cir- 
cles in Fig. |4] show estimates of Z for U/t = 6, which 
were obtained by fitting the y — 3toE(w„) data for the 
lowest three Matsubara frequencies x — uj n to a func- 
tion of the form x — Ay + By 2 . While a careful ex- 
amination of the behavior of Z near U = U C 2 may be 
desirable, these uncertainties do not affect the analysis 
in this paper. Using the approximation in Eq. (fT0|) . 
the (it = 400 data yield the specific heat coefficients 
(p/(U/2),jt) = (0.7, 38), (0.725, 45), (0.75, 59) for U/t = 
6 and (0.6, 33), (0.625, 43), (0.65, 61) for U/t = 6.5. 

Figure [5] shows the thermodynamic potential differ- 
ences between metallic and insulating solutions as a func- 
tion of chemical potential for U/t = 6 (upper panel) and 
U/t = 6.5 (lower panel). Because the thermodynamic 
potential differences are very tiny, taking proper account 
of the entropy of the metallic state is important. The 
point where the curve crosses zero yields the location of 
the first order transition, which is indicated by the verti- 
cal lines in Fig. [3] The phase transition occurs near the 
spinodal point and shifts with temperature in a similar 
way as the spinodal point. 

As a consistency check, we show the value of the dop- 
ing obtained from the thermodynamic potential using the 
formula 



2n v 



dSl n 



(12) 



as triangles in Fig. [3J These results, based on an approx- 
imation of the derivative by the slopes of the solid lines 
in Fig. [51 agree within 10% with the measured dopings. 
The thermodynamic potential curves at larger \fx — fi\\ 
show similar slopes for (it — 100, 200 and 400, and thus 
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FIG. 5: Solid points: thermodynamic potential difference be- 
tween the metallic and insulating solutions for the indicated 
temperatures, U = 6t (upper panel) and U — 6.5i (lower 
panel). Open symbols: "thermodynamic potential" computed 
by neglecting the 1 /2-yT 2 term in Eq. (J9]l . 



yield similar dopings. So, within the expected precision 
our thermodynamic potential analysis yields consistent 
results. 



B. Second order transition at T = 

To address the nature of the transition at T = we 
must first extrapolate the measured densities to the in- 



TABLE I: Location of the T — second order phase transition 
(fJ.ro), compressibility per spin dn/dfi, and coefficient of the 
quadratic term B(ji — /i C 2) 2 for different values of U/t. 
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FIG. 6: Close-up view of the small doping results for U > U C 2, 
showing the essentially linear onset of doping, which becomes 
more pronounced as one moves away from the critical point. 
The lines show parabolic fits to the data points which were 
extrapolated to T — > (U/t — 6, 6.5) or can be considered 
indistinguishable from that limit (U/t = 8). We assume that 
these curves correspond to n(T = 0, fi). 
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FIG. 7: Energy difference between the metallic and insulating 
solutions at T = 0, extrapolated from the data for fit = 200 
and fit = 400, respectively. The lines show the result obtained 
from the n(T = 0,/i), assuming a second order transition. 



sulating density n = 0.5. In the range of chemical po- 
tentials for which a metallic state is stable for both fi — 
200/t and (3 = 400/i we extrapolate the density to T = 
using the Fermi liquid relation n(T) = n(T = 0) + aT 2 
and fitting n(T — 0) and a. FigureOshows as solid points 
the result of the extrapolation to T = and as open sym- 
bols the computed density at our lowest temperature. 
One sees that in the density range (0.5 — n > 0.01) where 
more than one temperature is available, the — 400/t 
data are essentially converged to the T = value. The 
roughly linear dependence of doping on (/i — (Mi) in the 
region 0.5 — n < 0.03 — 0.05 (depending on U) implies a 
constant compressibility. Figure [HI shows that the density 
can be fitted very well over the entire measurement range 
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FIG. 8: Open circles: critical chemical potential fi C 2 deter- 
mined by the extrapolation of the density to temperature 
T = 0. Crosses: position of the spinodal point \iciifit = 400) 
(where the metallic solution ceases to exist) as a function of 
U. The stars indicate the size of the spectroscopic gap at 
half-filling, A, as determined by analytic continuation of in- 
sulating Green functions for fit = 40, and the solid line gives 
a rough estimate for this gap, which assumes a semi-circular 
density of states. 

to the function n(T = 0, fi) = A(fi — /i C 2) + B(fi — Hci) 2 ■ 
Performing the fit, we find the parameters listed in Tab.[U 
To verify the consistency of our analysis and determine 
the order of the transition, we show as lines in Fig. [7] the 
thermodynamic potential curves obtained by use of our 
fits to n(T = 0,(1) and the T = version of Eq. (fTT|) . 
with fi c set equal to the value /x C 2 at which n = 0.5 and 
the integration constant ri met (/u c ) set equal to E- msl 

E m M - £ins(^c2) -2 f dn'n(T = 0, //). (13) 

These curves are based on the assumption that at T = 
the energies of metallic and insulating states coincide 
only at the chemical potential /x = p, c i at which n = 0.5, 
and that the metal-insulator transition is hence continu- 
ous. 

On the other hand, we can extrapolate the measured 
thermodynamic potentials f£ me t(/i,T) at f3 = 400/t and 
/3 = 200/t to T = using the estimated 7 values and 
Eq. fJ5]) . These results are shown in Fig. [7] as solid and 
open points, respectively. The close agreement between 
the two estimates for the energy difference verifies the 
analysis and shows that at T = (unlike at T > 0) the 
doping does not jump discontinuously as the chemical 
potential is increased into the metallic region. 

The critical chemical potentials /i C 2 = fJ- C 2(T = 0) 
(measured from the half filling value /ii) are shown as 
open circles in Fig. [5] while the positions of the spin- 
odal points obtained from our solution at fi = 400/i are 
shown as crosses connected by the solid line. One sees 
that the zero temperature extrapolation is important for 
elucidating the behavior in the region near U C 2- We re- 
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mark that our T — > extrapolations agree well with the 
zero temperature results for /i c2 presented very recently 
in Ref . [Tfl 

It is interesting to compare the chemical potentials /i C 2 
with estimates of the gap at half filling. The stars in 
Fig. [8] indicate the size of the spectroscopic gap A ob- 
tained by analytically continuing Monte Carlo data for 
the higher temperature (3 — 40/t.— We see that the 
doping-induced states occur much before the chemical 
potential reaches the edge of the band; thus doping in- 
duces "in-gap" states. However, a glance at Fig. [3] shows 
that by the time the density is increased beyond a few 
percent, the chemical potential is inside the Hubbard 
bands. The in-gap nature of the states is therefore rel- 
evant only at extremely low dopings of a few percent or 
less. 

A scaling behavior is evident for interaction strengths 
near U C 2- In particular both the compressibility and the 
critical chemical potential vanish as U but the 

ratio remains roughly constant (see Tab. (Tp . A simple 
scaling analysis would suggest that 

n(/x)- 0-5 = n( f i- f iif$( ) (14) 

\Mc2 - Mi / 

with $(y) a function tending to a constant as y — > oo 
and vanishing at y — 1, but with a non- vanishing first 
derivative. The quadratic dependence of n on /i at larger 
chemical potentials and the rough scaling of <9n/cV| Mc2 
and (n C 2 — Mi) suggest an exponent x — 2. We do not 
have sufficient accuracy to determine precisely the scaling 
function and the behavior of fi C 2- These depend crucially 
on the value of U C 2, which we have not determined with 
precision. If the value U C 2 — 5.8t quoted in Ref. [H is 
used, our data are consistent with the relation 

( U — U \ 1 ^ 2 
Mc2 - fJ,i = Op — — . (15) 



IV. CONCLUSIONS 

We have studied the doping dependent Mott transi- 
tion in the one-band Hubbard model, using single site 
DMFT and a powerful diagrammatic QMC impurity 
solver which allows access to low temperatures even at 
strong interactions. A detailed quantitative understand- 
ing of the doping driven metal-insulator transition could 
be obtained. By computing the temperature dependence 
of the energy and thermodynamic potential we were able 
to perform a convincing extrapolation to T = and show 
that while the metal-insulator transition at T > is first 
order (with a jump in density), it becomes continuous 
at T = 0. At the critical chemical potential fi C 2(U), the 
density per spin, n, smoothly approaches the Mott insu- 
lating value 0.5. Our data are consistent with the scaling 
assumption that fi C 2(U) goes smoothly to the half filled 
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FIG. 9: dn/dfi normalized to the non- interacting value 
(dn/d(i)o ~ 1/nt, as a function of density per spin, 0.5 — n. 
For U close to U C 2, the compressibility is substantially reduced 
relative to the non-interacting system. 



band value [i\ = U/2 as U — > U^. Our results are in sub- 
stantial agreement with a very recently published density 
matrix renormalization group study of the same model.— 
This study determined (J, C 2(U) as the boundary of a coex- 
istence region, without making a statement on the order 
of the transition, finding for example (\i c i — /J-i){U — 
6t) w 0.5t and (\i c2 - (ii){U = 6.5t) « 0.9i, in good 
agreement with the estimates presented in Tab. |U 

We determined the behavior of the electronic com- 
pressibility dn/dfj, as a function of U and doping, finding 
that it vanishes at U — U C 2, [i = \i\ and grows roughly 
linearly with distance in U and fi from this critical point. 
The vanishing of the compressibility at the T > criti- 
cal end point of the Mott transition has been extensively 
discussed in the literature . 18 ! 19 In a series of interesting 
publications, Imada and co-workers have argued that at 
the density-driven T = metal-insulator transition the 
compressibility dn/d[i should vanish^ in contrast to our 
finding that the quantity is non- vanishing for U > U C 2- 
The conclusions of Imada and co-workers are based on 
hyperscaling, which is unlikely to apply in the limit of 
spatial dimensionality d —>■ oo in which the DMFT ap- 
proach is exact. Further consideration of this issue in 
finite dimensionality is an important open problem. 

The values we obtain for the electronic compressibil- 
ity are interesting. Figure [5] shows the compressibility, 
normalized to the non-interacting value of approximately 
1 /7rt, as a function of doping. These curves were obtained 
from the fitting functions for n(T = 0, fj,) and show that 
the suppression of the compressibility at ^ = Ml an d 
U = U C 2 persists over a wide interaction and doping 
range. This suppression has two experimental conse- 
quences: first, the square of the inverse Thomas-Fermi 
screening length 

2 47re 2 <9n 
<?tf = —^—q^ ( 16 ) 
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should be strongly reduced near the metal-insulator tran- 
sition, possibly leading to unusually weak screening of 
charged impurities. However, a simple estimate for 
high-T c materials gives dn/d/i w 1/eV. A lattice con- 
stant of 4A and a background dielectric constant of 10 

would then imply <7^p|band , so even the largest 

renormalization shown in Fig. [9] would only lead to a 
<7tf ~ 1A . Thus screening is always expected to be 
at the scale of a lattice constant. Nevertheless the effect 
might be observable in scanning tunneling microscopy. 
Possibly more easily observable would be a doping depen- 
dence of the sound velocity via the Bohm-Staver relation 
c 2 ~ dn/dfj,^ 

Finally, looking towards the future, we suggest that 
using the techniques presented here and in Ref. Il2| to 



reexamine the metal-insulator transition in the context 
of cluster dynamical mean field theories is an urgent open 
problem. 
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